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Search for progenitors of supernovae type la with SPY 
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Abstract. We have started a large survey for double degenerate (DD) binaries as potential progenitors of type la supernovae 
with the UVES spectrograph at the ESO VLT (ESO SN la Progenitor surveY - SPY). About 400 white dwarfs were checked 
for radial velocity variations during the first 15 months of this project, twice the number of white dwarfs investigated during 
the last 20 years. We give an overview of the SPY project and present first results Fifty four new DDs have been discovered, 
seven of them double lined (only 18 and 6 objects of these groups were known before, respectively). The final sample is 
expected to contain 150 to 200 DDs. Eight new pre-cataclysmic binaries were also detected. SPY is the first DD survey 
which encompasses also non-DA white dwarfs. SPY produces an immense, unique sample of very high resolution white 
dwarf spectra, which provides a lot of spin-off opportunities. We describe our projects to exploit the SPY sample for the 
determination of basic parameters, kinematics, and rotational vel ocities of white dwarfs. A catalogue with a first subset of 
our white dwarf data has already been published by Koester et al. (2001). 



1. Introduction 

Supernovae of type la (SN la) play an outstanding role for our 
understanding of galactic evolution and the determination of 
the extragalactic distance scale. However, the nature of their 
progenitors is still unknown (e.g. Livio 2000 1. There is gen- 
eral consensus that the event is due to the thermonuclear ex- 
plosion of a white dwarf when a critical mass (very likely the 
Chandrasekhar limit, 1.4A/©) is reached, but the nature of the 
progenitor system remains unclear. It must be a binary, with 
matter being transfered to the white dwarf from a compan- 
ion until the critical mass is reached. However, two options 
exist for the nature of the companion: either another white 
dwarf in th e so-c alled double degenerate (DD) scenario (Iben 
& Tutukov 1984), or a red giant/subgiant in the so-called sin- 



gle degenerate (SD) scenario (Whelan & Iben 1973 ), with the 
system possib ly appearing as a symbiotic binary (Munari & 
Renzini 1992 ) or a supersoft X-ray source (van den Heuvel et 



al. 1992) 



* Based on data obtained at the Paranal Observatory of the Euro- 
pean Southern Observatory for programs 165.H-0588, 167.D-0407, 
and 266.D-5658 

** Based on observations collected at the German-Spanish Astro- 
nomical Center, Calar Alto, operated by the Max-Planck-Institut fiir 
Astronomie Heidelberg jointly with the Spanish National Commis- 
sion for Astronomy 



The solution of the SD vs. DD dilemma is of great im- 
portance for assessing the role of SNe la in a variety of astro- 
physical situations, and - perhaps even more importantly - 
their eff ective ness as accurate cosmological probes (cf. Lei- 
bundgut 2000 ). As is widely known, current studies of high-z 
SNe la favo ur an accele rating , A-dominated universe (Riess 
et al. 1998 ; Leibundgut 2001 ). Yet, the result depends on the 
assumption that the light-curves of high-z SNe la are intrin- 
sically identical to those of the local SNe la. If instead high-z 
SNe la were just w 0.3 mag dimmer at maximum, then a 
A = universe would be preferred (Riess et al. 1998; Lei- 
bundgut 2000). Knowing whether the progenitors are SD or 
DD systems would constrain the possible explosion models 
and improve the use of SNe la as cosmological probes. 

Several systematic radial velocity (RV) searches for DDs 
have been undertaken starting in the mid 1980's (e.g. Robin- 
son & Shafter |1987| Bragaglia et al. |1990| Saffer, Livio & 



Yungelson 1998; Maxted & Marsh 1999). Before 2001, com 
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bining all the surveys, w 200 white dwarfs have been checked 
for RV variations with sufficien t accu racy yielding 18 DD s 
with periods P < 6 d .3 (Marsh |2000t Maxted et al. |2000a| ). 
None of the 18 systems seems massive enough to qualify as a 
SN la precursor. This is not surprising, as theoretical simula- 
tions suggest that only a few percent of all DDs are po tential 
SN la prog enitor s (Iben, Tutukov & Yungelson |1997| ; Nele- 
mans et al. 2001 ). 



In order to perform a definitive test of the DD scenario 
we have embarked on a large spectroscopic survey of 1500 
white dwarfs using the UVES spectrograph at the UT2 tele- 
scope (Kueyen) of ESO VLT to search for RV variable white 
dwarfs (ESO SN la Progenitor surveY - SPY). Details of 
our observations are given in Sect. ^. SPY will overcome the 
main limitation of all efforts so far to detect DDs that are 
plausible SN la precursors: the samples of surveyed objects 
were too small. 

The SPY program has now just reached mid term, but 
already dramatically increased the number of white dwarfs 
checked for RV variations and the number of known DDs. 
Our sample includes many short period binaries, several with 
masses closer to the Chandrasekhar limit than any system 
known before. An overview of our DD results is given in 
Sect.|. 

SPY produces a unique sample of high-quality, high- 
resolution spectra of white dwarfs. This is a large, homoge- 
neous database of enormous value for many other areas of re- 
search. Spin-off opportunities in general and projects started 
by us to exploit this treasure chest are discussed in Sect. || 

2. Observations and data reduction 

Targets for SPY are drawn from five sources: the white dwarf 
catalog of McCook & Sio n (|1999| ), the Hamburg ESQ Sur- 
vey (HES; Wisotzki et al. |2000| ; Chris tlieb e t al. |200l| ), the 
Hamburg Quasar Survey (Hagen et al. 



1995 



Homeier et al. 

(11998]), the M ontrea l-Cambridge-Tololo survey (MCT; Lam- 
ontagne et al. 200C), and the Edinburgh-Cape survey (EC; 
Kilkenny et al. 1997). Our selection criteria were spectro- 



scopic confirmation as white dwarf (at least from objective 
prism spectra) and B < 16.5. 

Spectra were taken with the UV- Visual Echelle Spectro- 
graph (UVES) of the UT2 telescope (Kueyen) of the ESO 
VLT. UVES is a high resolution Echelle spectrograph, which 
can reach a resolution of 1 10 ,000 i n the red region with a 
narrow slit (cf. Dekker et al. 2000 for a description of the 



instrument). Our instrument setup (Dichroic 1, central wave- 
lengths 3900 A and 5640 A) uses UVES in a dichroic mode 
with a 2048 x 4096 EEV CCD windowed to 2048 x 3000 
in the blue arm, and two CCDs, a 2048 x 4096 EEV and a 
2048 x 4096 MIT-LL, in the red arm. Nearly complete spec- 
tral coverage from 3200 A to 6650 A with only two «80A 
wide gaps at 4580 A and 5640 A is achieved. In the standard 
setting used for our observations UVES is operated with an 
8" decker in the blue arm and an 11" decker in the red arm. 

Our program is implemented as a service mode pro- 
gram. It takes advantage of those observing conditions, which 
are not usable by most other programs (moon, bad seeing, 



clouds) and keeps the VLT busy when other programs are 
not feasible. A wide slit (2.1") is used to minimise slit losses 
and a 2 x 2 binning is applied to the CCDs to reduce read 
out noise. Our wide slit reduces the spectral resolution to 
R = 18 500 (0.36 A at Ha) or better, if seeing disks were 
smaller than the slit width. Depending on the brightness of 
the objects, exposure times of 5 min or lOmin were chosen. 
The S/N per binned pixel (0.03 A) of the extracted spectrum 
is usually 15 or higher. Due to the nature of the project, two 
spectra at different, "random" epochs separated by at least 
one day are observed. 

Although our program is carried out during periods of less 
favourable observing conditions, the seeing is often smaller 
than the selected slit width of 2.1". This can, in principle, 
cause wavelength shifts, if the star is not placed in the center 
of the slit. However, since according to the standard observ- 
ing procedure the star is first centered on a narrow slit before 
the exposure with the broader slit is started, it can be expected 
that the star is usually well centered in the slit. We note that 
possible wavelength shifts can be corrected with the telluric 
absorption features present in the red region of our spectra. 

Data reduction: ESO provides a data reduction pipeline 
for UVES, based on MIDAS procedures. The raw data are 
bias and interorder background subtracted and extracted with 
an optimum extraction algorithm. The orders are flatfielded, 
rebinned, and merged. Wavelength calibration is performed 
by means of ThAr calibration lamp exposures. The quality of 
the reduced spectra is in most cases very good; especially the 
removal of the interorder sensitivity variation and merging of 
the orders works very well. Sometimes the reduction pipeline 
produces artifacts of varying strength, e.g. a quasiperiodic 
pattern in the red region similar in appearance to a fringing 
pattern. In a few cases either the blue or the red part of the 
spectrum has extremely strong artifacts of unknown orign. A 
"UVES" context for MIDAS is available, which provides re- 
duction routines adapted for the UVES spectrograph and can 
be used to re-reduce spectra. As an example of the quality 
achievable with the pipeline reduction one spectrum without 
extraction problems is shown in Fig. [|. 

Remaining large scale variations of the spectral response 
function are removed by utilizing nearly featureless spectra of 
DC white dwarfs or of sdO stars with narrow spectral lines, 
which are normalised and used to derive an approximate re- 
sponse function. This function was applied to the other spec- 
tra. Since the sampling of the spectra is much higher than 
needed for many purposes (besides measuring RVs), we pro- 
duced a version of the spectra, which were rebinned to 0. 1 A 
stepsize and smoothed with a Gaussian of 0.2 A FWHM. This 
produces only a slight degradation of the resolution while 
considerably improving the signal-to-noise level. 

Radial velocity measurement: Since SPY produces a 
large number of spectra, which have to be checked for RV 
variations, a fast and reliable algorithm to measure RV shifts 
is necessary. One often applied standard technique is cross 
correlation of observed spectra. The RV shift is measured 
from the peak of the correlation function. However, it turned 



R. Napiwotzki et al.: Search for progenitors of SNIa with SPY 



3 




6600 



Fig. 1. Pipeline reduced UVES spectrum of the DA white dwarf HE 0403— 4129. The panels correspond to the blue channel 
and to both parts of the red channel covered by different CCDs. Artifacts are marked by asterisks. The spectrum was slightly 
smoothed and rebinned as described in the text. 



out that this method is not well suited for the determination 
of RVs from the narrow cores of the broad Balmer lines of 
DA white dwarfs. This can be understood as follows. If we 
limit the correlation to, say, the inner 10 A of the line to avoid 
spurious signals caused by normalization problems etc., the 
boundaries of the integration interval are located in the line 
wings. Contributions from both interval boundaries (which 
generally do not cancel out) result in a large extra signal in 
addition to the real correlation signal, which makes the deter- 
mination of the correlation peak very difficult or even impos- 
sible. 

We overcome this problem by applying a correlation 
based on a x 2 test instead of the usual cross correlation. Our 
procedure starts with a noise estimate from the rms scatter 
computed within a moving box. Cosmics etc. are flagged out 
with an clipping algorithm. Since it is important for comput- 
ing a realistic \ 2 that the spectra have the same flux level, we 



scale the second spectrum to the first one using a low order 
polynomial to account for possible differences of the contin- 
uum slope. Afterwards \ 2 1S calculated as a function of RV 
difference, by shifting the second spectrum in small RV steps 
(0.1km s _1 ) and calculate x 2 by direct numerical integra- 
tion (without rebinning). Since the calculation of \ 2 f° r two 
noisy spectra results in some noise in the resulting x 2 curve, 
finally a slight smoothing (Gaussian of FWHM 1km s _1 ) 
is applied. The RV shift is evaluated from the minimum % 2 . 
Error margins can be estimated by determining the intervals 
corres ponding to a certain A\ 2 as described in Press et al. 
( |1992|) . 

From a formal point of view x 2 tests were invented to 
test and find best fitting models. Nevertheless, our method 
works and produces reliable results with the observed spec- 
tra as well. We performed several tests: Two observed spec- 
tra were correlated with a model spectrum and the RV shifts 
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Table 1. Fraction of RV variable stars in the current SPY sam- 
ple for different spectral classes. WD+dM denotes systems 
for which a previously unknown cool companion is evident 
from the red spectra. 



Spectral type 


total 


RV variable 


detection rate 


All 


387 


62 


16% 


DA 


310 


53 


17% 


non-DA (DB,DO,DZ) 


53 


1 


2% 


WD+dM 


24 


8 


33% 



added, which yielded virtually the same result as the direct 
X 2 correlation of both observed spectra. We also made exper- 
iments with simulated spectra, produced from model spectra 
by adding artificial noise and known RV shifts. These RV dif- 
ferences could always be reproduced from the \ 2 correlation 
within the estimated error limits. One great advantage of our 
procedure is its flexibility and that it can easily be applied to 
measure RV shifts in stars of different spectral types (Balmer 
lines of DA white dwarfs, He I lines of DBs, and metal lines 
of sdB stars like HE 1047—0436 discussed below). We rou- 
tinely measure RVs with an accuracy of « 5 km s - 1 or better, 
therefore running only a very small risk of missing a merger 
precursor, which have orbital velocities of 150km s _1 or 
higher. 



3. Results on double degenerates 

We have analysed spectra of 387 white dwarfs taken during 
the first 15 months of the SPY project and detected 54 new 
DDs, 7 are double-lined systems (only 6 were known before). 
The great advantage of double-lined binaries is that they pro- 
vide us with a well determined total mass. Since it is likely 
that the SPY sample contains even more double-lined sys- 
tems (with a faint secondary), we will check follow-up obser- 
vations of apparently single-lined systems for the signature of 
the secondary. Our sample includes many shor t period bina- 
ries (some examples are discussed in Sect. 3.1), several with 
masses closer to the Chandrasekhar limit than any system 
known before. In addition, we detected 8 RV variable systems 
with a cool main sequence companion (pre-cataclysmic vari- 
ables; pre-CVs). Some examples of single-lined and double- 
lined DDs are shown in Fig. |2|. A more detailed inventory of 
the SPY sample is given in Table [j] Our observations have 
already increased the DD sample by a factor of four. After 
completion, a final sample of 150 to 200 DDs is expected. 

Follow-up observations of this sample are mandatory to 
exploit its full potential. Periods and white dwarf parameters 
must be determined to find potential SN la progenitors among 
the candidates. Good statistics of a large DD sample will also 
set stringent constraints on the evolution of close binaries, 
which will dramatically improve our understanding of this 
phase of stellar evolution. During our follow-up observations 
we have detected a very promising potential SN la precursor 
candidate. However, some additional spectra are necessary to 
verify our RV curve solution. Results will be reported else- 
where. 




MA) 

Fig. 2. Three double-lined and one single lined RV variable 
DDs from our VLT survey. The dotted line marks the rest 
wavelength of Ha. The spectra are rebinned and smoothed to 
a resolution of «0.5 A. 



Although important information like the periods, which 
can only be derived from follow-up observations, are 
presently lacking for most of the stars, the large sample size 
already allows us to draw some conclusions. (Note that fun- 
damental white dwarf parameters like masses are known from 
the spectral analysis described below). One interesting as- 
pect concerns white dwarfs of non-DA classes (basically the 
helium-rich spectral types DB, DO, and DZ, in contrast to 
the hydrogen-rich DAs). SPY is the first RV survey which 
performs a systematic investigation of both classes of white 
dwarfs: DAs and non-DAs. Previous surveys were restricted 
to DA white dwarfs, because the sharp NLTE core of Ha pro- 
vides a very accurate RV determination. This feature is not 
present in non-DA (DB, DO) spectra, but the use of several 
helium lines enables us to reach a similar accuracy. Thus the 
low number of DDs found for the non-DA spectral types is 
statistically significant (Table [l]) and indicates that important 
evolutionary channels for the formation of DDs produce pref- 
erentially hydrogen-rich white dwarfs of the DA variety. 

3.1. Highlights of our follow-up results 

Follow-up observations are necessary to determine the sys- 
tem parameters of the DDs. We concentrated on candidates 
with high RV variations, indicating short periods, because 
the probability to find potential SN la candidates among these 
systems is highest. However, let us note that probably some 
of the "small ARV" DDs could be short period systems (pos- 
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sibly even SNIa progenitors) with low inclination angles 
and/or unfavourable phase differences of the SPY observa- 
tions. Follow-up observations have been carried out at the 
VLT as well as with the 3.5 m telescope of the Calar Alto 
observatory /Spain. 

In the following paragraphs we present results for two of 
our most interesting binaries: the subdwarf B + white dwarf 
system HE 1047-0436 (Napiwotzki et al. |2001a| ) and the 
do uble-lin ed DA+DA binary HE 1414-0848 (Napiwotzki et 
al. 2001b) with a mass only 10% below the Chandrasekhar 
limit. 



HE 1047-0436: HE 1047-0436 was discovered by the 
Hamburg ESO Survey (HES) as a potential hot white dwarf 
and, therefore, included in our survey. The UVES spectra 
(Fig. ||) showed that it is in fact a subluminous B star (sdB) 
with a rather large RV shift of 160 km s _1 , which made the 
star a prime target for further study. SdB stars are pre-white 
dwarfs of low mass (w 0.5M©) still burning helium in the 
core, which will evolve directly to the white dwarf stage omit- 
ting a second red giant phase. 

Recently, sdBs with white dwarf components have been 
propose d as potential SNela progenitors by Maxted et al. 
( |2000b| ) who discovered that KPD 1930+2752 is a sdB+WD 
system. Its total mass exceeds the Chandrasekhar mass and 
the components will merge within a Hubble time, which 
makes KPD 1930+2752 a SNIa progenitor candidate (al- 
though this inter pretation has been questioned recently, 
Ergma et al. |200l| ). 

The high RV shift of HE 1047-0436 prompte d us t o 
perform follow-up observations (Napiwotzki et al. 2001a). 
The orbital period of P = 29 h 7 m 5 s , a semi-amplitude of 
94 km s _1 , and a minimum mass of the invisible companion 
of 0.44 M are derived from the analysis of the RV curve 
(Figs. Q and ||). We use an upper limit on the projected ro- 
tational velocity of the sdB star to constrain the system in- 
clination and the companion mass to M > 0.71 M , bring- 
ing the total mass of the system closer to the Chandrasekhar 
limit. However, the system will merge due to loss of angular 
momentum via gravitational wave radiation only after sev- 
eral Hubble times. Atmospheric parameters (T e ff = 30200 K, 
logg = 5.7) and metal abundances are also derived (cf. 
Fig. ||). The res ulting values are typical for sdB stars (cf. e.g. 
Edelmann et al. 2002). 



HE 1414-0848: HE 1414-0848 is a double-lined 
DA+DA binary (Fig. ||). The orbital period of 
P = 12 h 25 m 44 s and semi-amplitudes of 127 km s" 1 
and 96kms~ 1 are derive d for the individual components 
(Napiwotzki et al. 2001b ). RV curves for both components 
are displayed in Fig. |7j. The ratio of velocity amplitudes 
is directly related to the mass ratio of both components. 
Additional information comes from the mass dependent 
gravitational redshift z = AA/A = GMR^ x cr 2 , which 
for a given mass-radius relation can be computed as a 
function of white dwarf mass. The difference in gravitational 
redshift corresponds to the apparent difference of "systemic 
velocities" of both components, as derived from the RV 
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Fig. 3. Fit of the hydrogen and helium lines of the RV cor- 
rected, coadded UVES spectra of HE 1047-0436. 
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Fig. 4. Power spectrum of the HE 1047-0436 measurements 
used for the periodogram analysis. 



curves. Only one set of individual white dwarf masses fulfils 
the constraints given by both the amplitude ratio and redshift 
difference. We estimate the masses of the individual com- 
ponents with this method: 0.54AT Q and 0.72M Q . Another 
estimate of the white dwarf parameters is available from 
a model atmosphere analysis of the combined spectrum. 
The formal result indicates an average mass of O.68M 
consistent with the results derived from the analysis of the 
RV curve. The total mass of the HE 1414-0848 system is 
1.26M Q , only 10% below the Chandrasekhar limit! The 
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Fig. 5. Measured RVs as a function of orbital phase and fitted 
sine curve for HE 1047-0436. 
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Fig. 6. Ha spectra of HE 1414-0848 covering 5 hours during 
one night together with a fit of the line cores. The numbers 
indicate the Julian date of the exposures and the orbital phase 



system will merge due to loss of angular momentum via 
gravitational wave radiation after two Hubble times. There- 
fore HE 1414—0848 does not qualify as a SNIa progenitor, 
but it is the most massive short period DD known today. 




0.5 
Phase 



Fig. 7. Measured RVs as a function of orbital phase and fit- 
ted sine curves for HE 1414-0848. Filled circles/open rect- 
angles indicate the less/more massive component. Note the 
difference of the "systemic velocities" 70 between both com- 
ponents caused by gravitational redshift. 



4. Spin-off results 

SPY produces an immense, unique sample of very high res- 
olution white dwarf spectra. This database will have a large 
impact on many fields of white dwarf science. It will allow 
us for the first time to tackle many longstanding questions on 
a firm statistical basis. Among those are the mass distribution 
of white dwarfs, kinematical properties of the white dwarf 
population, surface compositions, luminosity function, rota- 
tional velocities, and detection of weak magnetic fields. We 
have therefore decided to make these results available to the 
community in the form of a catalog with some preliminary 
interpretation as soon as feasible, not waiting for the com- 
pletion of the SPY project. A first part of thi s cata log was 
published in the recent paper of Koester et al. (2001 ), cover- 
ing observations of about 200 white dwarfs of spectral types 
DA and DB. For all spin-off opportunities mentioned above 
the statistics will be dramatically improved by the final white 
dwarf spectra database. We are exploiting the SPY sample 
for two spin-off projects, which take advantage of the high 
spectral resolution: the kinematics of white dwarfs and their 
rotational velocities. 



Basic parameters and mass distribution of white dwarfs: 

The response of the UVES spectrograph is sufficiently stable 
and well behaved to allow the determination of temperatures, 
gravities, and masses for the programme white dwarfs. Thus 
we carried out a model atmosphere analysis and p ublished 
the results in our catalog paper (Koester et al. 2001). Agree- 
ment of our results with literature values is quite good for 
the DA sample, thus confirming that the UVES spectra can 
be used for accurate and reliable parameter determination of 
white dwarfs. This was not necessarily expected, since the in- 
dividual Balmer lines in the DA span several echelle orders. 
Note that our catalog paper included the second ever pub- 
lished analysis of a sample of helium-rich DB white dwarfs. 
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Fig. 8. Meridional projection of the Galactic orbits of 
WD 0013-241 (top) and WD 0016-258 (bottom) computed 
over lOGyr. WD 0016—258 is a typical thin disk star with an 
orbit of small eccentricity reaching only small heights above 
the Galactic plane. The orbit of WD 0013-241 on the other 
hand has high eccentricity and reaches larger distances from 
the Galactic plane, typical for a thick disk star. 



Our parameter determinations allow the selection of spe- 
cial types of white dwarfs for further follow-up investigations 
(by us as well as by other members of the community), and 
are a prerequisite for the determination of space velocities 
(see next paragraph). 




xrAi 

Fig. 9. Observed core of the Ha line compared to synthetic 
spectra computed for different rotational velocities vsmi. 
This spectrum allows to put an upper limit of 15 km s _1 on 



Examples of thin and thick disk orbits are shown in Fig. ||| 
So far 47 white dwarfs have been analysed. This will be the 
first kinematic study of a large sample of white dwarfs, which 
is based on the complete set of space motion in three dimen- 
sions. 

Rotational velocities: High resolution spectra of the Ha 
line are very useful for the determination of rotational ve- 
locities (v sin i) of white dwarfs, which puts constraints on 
the evolution of angular moment um du ring stellar evo lution . 
Two recent studies (Heber et al. 1997, Koester et al. 1998) 



increased the number of determinations just to 54. It turned 
out that white dwarfs are generally very slow rotators with a 
few, puzzling exceptions. 

Many spectra taken for the SPY project allow the determi- 
nation of rotational velocities with better accuracy than previ- 
ous studies. We started the systematic analysis of the present 
sample (one example is displayed in Fig. ^). The final sample 
will dramatically increase the number of white dwarfs with 
known rotational velocities and allow for the first time to ex- 
plore trends of v sin i with, e.g., cooling age. 



Kinematical properties of white dwarfs: The large SPY 
sample of white dwarfs can be used to study the kinematics 
of these stars in our Galaxy, especially with respect to their 
belonging to the thin or thick disk population. The aim is to 
verify recent claims that a substantial fraction of the mass of 
the Galactic disk may be provided by thick disk white dwarfs, 
as it is concluded from the abundance analysis of a l ocal sa m- 
ple of F and G type main sequence stars (Fuhrmann 2001 ). 

We combined the measured RVs of the white dwarfs and 
determined proper motions from the DSS1/DSS2 and USNO 
data. Distances and corrections of the RV for gravitational 
redshift are computed from the fundamental parameters de- 
rived by Koester et al. (2001 ). Thus the space motion in three 



dimensions is known and Galactic orbits can be computed 



Chemical abundances: High-resolution spectra are very 
valuable for the determination of surface compositions. Most 
white dwarfs have very pure atmospheres containing only hy- 
drogen (DA) or helium (DO, DB). However, in the atmo- 
spheres of some hotter white dwarfs as well as in a frac- 
tion of cooler white dwarfs (types DZ, DAZ, DBZ) metal 
lines can be detected. Our understanding of this phenomenon 
is currently limited by the insufficient sample size of white 
dwarfs observed with high enough resolution. The largest 
published sample contains 38 white dwarfs observed with the 
HIRES spectrograph of the Keck I telescope (Zuckerman & 
Reid 1998). It is obvious that SPY can help to improve the 



from a code developed by Odenkirchen & Brosche (1992) 



statistics by a large factor. 

A few of the white dwarf candidates observed by SPY 
turned out to be subdwarf B stars (i.e. pre-white dwarfs with 
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Fig. 10. Chemical abundances in the atmosphere of the sub- 
dwarf B star HE 1047— 0436 compared to the solar values. 
Numbers indicate the number of lines measured for each ion- 
isation stage. 

lower surface gravities). These stars display puzzling abun- 
dance patterns likely caused by diffusion processes (cf. Edel- 
mann et al. 2002). As an example for the resulting chemical 
abundance pattern we present in Fig. [lO] our analysis of the 
binary HE 1047-0436 discussed above (Sect. O). 



5. Summary 

SPY has now reached about mid term and it has already 
tripled the number of white dwarfs checked for RV variablity 
(from 200 to 600) and quadrupled the number of known DDs 
(from 18 to 72) compared to the results achieved during the 
last 20 years. Our sam ple includes m any short period bina- 
ries (Napiwotzki et al. 2001a 2001b ), several with masses 
closer to the Chandrasekhar limit than any system known be- 
fore, greatly improving the statistics of DDs. We expect this 
survey to produce a sample of 150 to 200 DDs. 

This will allow us not only to find several of the long 
sought potential SN la precursors (if they are DDs), but will 
also provide a census of the final binary configurations, hence 
an important test for the theory of close binary star evolution 
after mass and angular momentum losses through winds and 
common envelope phases, which are very difficult to model. 
An empirical calibration provides the most promising ap- 
proach. A large sample of binary white dwarfs covering a 
wide range in parameter space is the most important ingredi- 
ent for this task. 

SPY produces a unique sample of white dwarf spectra 
with many spin-off opportunities, which will have a large im- 
pact on white dwarf science. We published a model atmo- 
sphere analysis of a first set of 200 white dwarfs (Koester et 
al. 2001 ). This allowed us to construct a mass distribution of a 



large sample of white dwarfs and made this sample known to 
the community for further follow-up investigations. We are 
exploiting the SPY sample for two spin-off projects, which 
take advantage of the high spectral resolution: we study the 
kinematics of white dwarfs and determine their rotational ve- 
locities. Members of the community interested in spin-off op- 



portunities are invited to participate in the exploitation of the 
SPY sample. 
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